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ABSTRACT 

Starting from the quasar sample of the Sloan Digital Sky Survey (SDSS) for which 
the CiV line is observed, we use an analysis scheme to derive the z-dependence of the 
maximum mass of active black holes, which overcomes the problems related to the 
Malmquist bias. The same procedure is applied to the low redshift sample of SDSS 
quasars for which H/3 measurements are available. Combining with the results from 
the previously studied Mgn sample, we find that the maximum mass of the quasar 
population increases as (x+z) 1 ' 64 ™' 04 in the redshift range 0.1 ^ z 4, which includes 
the epoch of maximum quasar activity. 

Key words: galaxies: evolution - galaxies: active - galaxies: nuclei - quasars: general. 



1 INTRODUCTION 

The dependence of the quasar black hole mass (Mbh) on 
redshift is of fundamental cosmological interest, since it is 
a direct probe of the supermassive black hole (BH) forma- 
. tion and evolution and of the triggering of the active galac- 
' tic nuclei (AGN) phenomenon. A full understanding of the 
cosmic evolution of quasar population would also give im- 
portant insights into the joint formation of BHs and their 
host galaxies. 

The basic procedure for the determination of the BH 
mass of AGN is through a direct application of the virial the- 
orem. In the hypothesis that the clouds responsible for broad 
lines emission move in Keplerian orbits (see e.g. Peterson & 
Wandel 2000 for supporting evidence), the mass depends on 
the clouds velocity, which is constrained by the line-width, 
and the cloud distance, that was shown to be strictly linked 
to the continuum intensity (e.g. Kaspi et al. 2000; Vester- 
gaard 2002; Kaspi et al. 2005; see also references in Section 
2). 

In order to explore the Mbh (z) dependence, one should 
take into account the Malmquist bias, related to the mini- 
mum detectable flux, which affects all the large samples of 
quasars and is apparent as an increase of the average ob- 
served luminosity with z. As Mbh is related to the nuclear 
luminosity, a Malmquist-type bias also affects the apparent 
trend of quasar BH masses (see e.g. Vestergaard et al. 2008 
and Kelly, Vestergaard & Fan 2009 for a detailed bias anal- 
ysis). 

In a previous paper (Labita et al. 2009, hereafter Pa- 
* E-mail: marzia.labita@gmail.com 



per I) we proposed a procedure to trace the z-dependence 
of Mbh overcoming the problems related to the Malmquist 
bias. We considered the ~ 50000 quasars from the Sloan 
Digital Sky Survey (SDSS) Data Release 5 (DR5, Schnei- 
der et al. 2007; see also Richards et al. 2002) for which Mgn 
line-width (FWHM) and 3000 A continuum luminosity were 
measured (Shen et al. 2008). Our scheme requires to fit the 
distribution of quasars in the luminosity-FWHM plane (see 
also Fine et al. 2008), using as free parameters the max- 
imum mass, the Eddington ratio limit and the minimum 
detectable luminosity. We obtained a dependence with red- 
shift of the type log(MBH X /M ) ~ 0.3z + 9 in the interval 
0.35 is £ ~ 2.25, which indicates that most massive BHs ex- 
hibit quasar activity at smaller cosmic times than less mas- 
sive ones (downsizing, e.g. Merloni 2004 and Shankar, Wein- 
berg & Miralda-Escude 2009 for a theoretical approach) . 

Here we apply the analysis described in Paper I on the 
SDSS quasar sample for which Civ based virial BH mass es- 
timates can be derived (Section l3.2p . This extends the study 
of the z-dependence of the mass of active black holes to 
2 ~ 4, including the redshift region where the quasar activ- 
ity is maximal (zp Ca k ~ 2 — 3, e.g. Richards et al. 2006a). To 
improve the coverage of the low redshift region, which is only 
marginally described by Mgn data, the adopted procedure 
is also applied to the SDSS quasars for which the virial BH 
mass can be derived from measurements of the H/3 spectral 
range (Section 13. 3|) . 

The dependence on redshift of quasar BH mass from 
the H/3, Mgn and Civ analysis is derived in Section [4] and 
the results are discussed in Section [S] To the best of our 
knowledge, this is the first time the Mbh{z) dependence 
of active black holes is derived beyond the peak of quasar 
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Table 1. Rcdshift coverage of the H/3, Mgn and Civ samples. 



Sample 


Objects 


■^min 




■2- max 


H/3 


11338 


0.08 


0.54 


0.91 


Mgn 


46649 


0.32 


1.27 


2.38 


Civ 


14070 


1.79 


2.42 


4.98 
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Figure 1. Distribution in redshift of the H/3 sample (red thin 
shading), of the Mgn sample (green medium shading) and of the 
Civ sample (blue thick shading). 

activity from a large sample of AGN, taking into account 
the problems related to the Malmquist bias. 

Throughout this paper, we adopt a concordant cosmol- 
ogy with Ho = 70 km s -1 Mpc -1 , Q. m = 0.3 and Q. A = 0.7. 
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Figure 2. The 8 panels show the Civ sample in the FWHM- 
luminosity plane at increasing redshift. Dotted, dashed and dash- 
dotted lines represent the loci of constant monochromatic lumi- 
nosity, constant mass and constant Eddington ratio respectively. 

Civ sample and the calibration of McLure & Dunlop (2004) 
for the H/3 sample: 

log .Mbh = 6 + log(a) + 21og(FWHM) + b log \L\ (1) 

where a = 4.57 and b — 0.53 with the Civ line-width and 
the 1350 A luminosity and a = 4.70 and b = 0.61 with 
the H/3 line-width and the 5100 A luminosity. Here Mbh is 
expressed in solar masses, FWHM in units of 1000 km/s and 
\L\ in units of 10 44 erg/s. To derive the Eddington ratio 
( £ t " ?1 ) from spectroscopic data, we adopt the bolometric 
corrections by Richards et al. (2006b). 



2 THE Civ AND H/3 SAMPLES 

Shen et al. (2008) studied ~ 80000 quasars included in the 
SDSS DR5 quasar catalogue. For ~ 60000 of them, they pro- 
vide rest-frame line widths and monochromatic luminosities 
at 5100 A, 3000 A and 1350 A from which it is possible 
to evaluate the BH mass through the virial theorem (see 
Shen et al. 2008 for details on calibrations and procedures 
of spectral measurements). 

In Paper I, we analysed the so-called Mgn sample 
(0.35 ~ 2 ~ 2.25), containing ~ 50000 quasars from Shen et 
al. (2008) with Mgn measures. Here we focus on the ClV 
sample, consisting of the ~ 15000 quasars from Shen et 
al. (2008) for which Civ line width and 1350 A monochro- 
matic flux are available. The analysis will be limited to ob- 
jects with z < 4, since statistics are extremely poor at higher 
redshift. We also consider the low redshift (0.1 is z ^ 0.9) 
SDSS DR5 quasars with readily available H/3 line width and 
5100 A flux (~ 10000 objects; hereafter, H/3 sample). All 
together, the H/3, Mgn and Civ samples cover ~ 60000 ob- 
jects and spread over a wide redshift range, from z w 0.1 to 
2 ?a 4 (see Table [Q and Fig. [TJ|. 

Following Shen et al. (2008), we assume the applicabil- 
ity of the virial theorem and adopt the calibration of Vester- 
gaard & Peterson (2006) to evaluate the BH mass for the 



3 THE MAXIMIZING MASS PROCEDURE 

3.1 Description of the basic procedure 

In Paper I we developed a procedure to derive the Mbh(z) 
trend overcoming the problems related to the Malmquist 
bias, for the Mgn sample of SDSS quasars. The objects 
were divided in 8 redshift bins with constant co-moving vol- 
ume. In each, we considered the objects distribution in the 
FWHM-luminosity plane. We aimed to construct a proba- 
bility density to describe the observed distribution in terms 
of a minimum luminosity (Z m in) due to the survey flux limit, 
a maximum mass (m ma x) and a maximum Eddington ra- 
tio (e max ; by definition, I = logXLx, m = log and 

e = log ^ bo1 ). These three cuts form a triangle, which de- 

^Edd ' ° 

scribes qualitatively well the shape of the quasar distribution 
in the FWHM-luminosity plane (see fig. 2 of Paper I) . Then, 
in each redshift bin, the assumed probability density (Eq. 4 
of Paper I) and the observed distribution of objects were dis- 
cretized in 600 boxes, and they were compared with a best- 
fitting procedure in order to determine the free parameters 
(^min, Wmax and e max and the widths of the corresponding 
distributions ai, a m and a e ) as a function of z. Errors on the 
best-fitting parameters were determined with Monte Carlo 
simulations (see Paper I for details on the procedure). 
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Figure 3. The 8 panels show the Civ sample in the FWHM— luminosity plane at increasing redshift: solid black contour plot (levels: 
30, 80, 140, 250, 400, 750 objects per box, see text) represents the discrete observed distribution of objects. Dotted red contour plot 
(same levels) shows the discrete distribution of a sample of objects simulated with the Monte Carlo method, adopting the assumed 
Pi fwhm(^, FWHM) probability density (Eq. 4 of Paper I) with the best-fitting parameters (see Table [2}. Dotted, dashed and dash- 
dotted lines represent Z m i n , m max and e m iv a x (fixed) respectively. In each panel, the shaded area shows the region removed from the fit. 
The rectangle in the first redshift bin represents the boxing adopted to discretize the distributions. 



Table 2. Best-fitting values of minimum luminosity, maximum mass and widths of the 
corresponding distributions, with errors and xtt f° r * ne Civ sample. In all the redshift bins, 
the number of degrees of freedom is v = 46 (50 data-points and 4 free parameters). 



Bin 


< z > 


#(Objects) 


^min 


07 


^max 


Cm 


xl 


1st 


1.94 


3484 


45.57±0.06 


0.19±0.03 


9.52±0.08 


0.42±0.05 


1.55 


2nd 


2.20 


1988 


45.74±0.07 


0.27±0.05 


9.50±0.10 


0.42±0.05 


2.44 


3rd 


2.46 


982 


45.86±0.10 


0.25±0.08 


9.61±0.07 


0.40±0.06 


2.22 


4th 


2.73 


325 


45.94±0.08 


0.23±0.05 


9.63±0.08 


0.44±0.05 


1.33 


5th 


3.02 


643 


45.95±0.08 


0.24±0.06 


9.68±0.08 


0.43±0.07 


2.10 


6th 


3.27 


605 


45.95±0.10 


0.29±0.06 


9.71±0.09 


0.48±0.08 


1.74 


7th 


3.58 


362 


46.06±0.09 


0.27±0.06 


9.72±0.11 


0.46±0.08 


1.00 


8th 


3.83 


236 


46.20±0.08 


0.23±0.06 


9.73±0.11 


0.45±0.07 


1.32 



Note. Data that come from a best-fitting procedure are displayed in boldface. 



3.2 Analysis of the Civ sample 

We aim to reproduce our analysis scheme on the Civ sam- 
ple. A key point is that the Civ emitting region is disc-like 
(e.g. McLure & Dunlop 2002; Labita et al. 2006, Decarli et 
al. 2008), as it is suggested by considerations on the broad 
lines shape and from the geometrical factor /, which links 
the observed FWHM of the line to the virial velocity (see 
McLure & Dunlop 2002 for definition and discussion on the 
geometrical factor). The emission line broadening is due to 



Doppler effect, and it is related only to the velocity compo- 
nent in the line of sight direction. If the Civ broad line region 
(BLR) is disc-like and it is observed face-on, the line width 
is null independently of the virial velocity. A consequence of 
this is that the data distribution in the FWHM-luminosity 
plane (Fig. [2]| is stretched towards lower values of the line 
width, with respect to the distribution that we would obtain 
if the BLR was isotropic. 

This implies that the observed distribution of objects 
in the FWHM-luminosity plane is not well described by the 



4 M. Labita, R. Decarli, A. Treves and R. Falomo 



assumed probability density (Eq. 4 of Paper 1), which is con- 
structed under the assumption that FWHM is proportional 
to the BLR virial velocity. In fact, if we try to reproduce 
the fit procedure described in Paper I to the Civ sample, we 
obtain that the \v values of the best fits are much bigger 
than those obtained in the analysis on the Mgn sample. In 
particular, the description results mediocre for low values of 
the FWHM, while the minimum luminosity and the maxi- 
mum mass sides of the triangle still describe adequately the 
data distribution, confirming that it appears stretched to- 
wards lower values of FWHM as expected. To rule out the 
possibility that this is due to bad spectral measurements, we 
checked the FWHM values provided by Shen et al. (2008) 
for a randomly selected subsample of quasars with narrow 
Civ lines, finding a very good agreement. 

In order to overcome the problem, we restrict our fit pro- 
cedure to the region of the FWHM-luminosity plane with 
FWHM>5000 km/s, which corresponds to those objects 
that are observed at high inclination angles (8 > 20° — 25° 
where 8 is the angle between the line of sight and the normal 
to the disc plane, Decarli et al. 2008). Of course there will 
still be some orientation effects in the FWHM>5000 km/s 
sample. However, the key point is that at these line-widths 
the FWHM distribution predicted under isotropy assump- 
tions is practically indistinguishable from that predicted in 
a disc-like BLR picture (see again Decarli et al. 2008, and 
in particular fig. 6 of their paper). 

Using this approach we lose the information on the max- 
imum Eddington ratio, since the parameters e max and a e 
are no more constrained as they influence the shape of the 
probability density mostly in the lower part of the FWHM- 
luminosity plane. We thus fix the values of these parameters 
to those obtained from the Mgn sample (e max = —0.34 and 
(j e = 0.22). This choice is supported by the fact that they 
appear to be independent of z in the redshift range studied 
through the Mgn data (0.35 < z < 2.25, see Paper I). 

We then apply the fit procedure described in Paper I 
on the high FWHM region of the FWHM-luminosity plane 
with 4 free parameters (Z m i n , o";, ?rimax and er m ) to be con- 
strained in all the redshift bins. The value of Z m i n is not con- 
strained by cosmology as in Paper I, because the SDSS DR5 
flux limit is not homogeneous. This effect is particularly pro- 
nounced at z ~ 3, because the selection algorithm basically 
targets ultraviolet excess quasars to i = 19.1 and higher red- 
shift (z > 3) quasars to i = 20.2 (see Richards et al. 2002 and 
http://www.sdss.org/dr5/algorithms/target.html for a 
detailed description of the SDSS target selection). 

The expected and the observed distribution of objects 
are discretized in boxes with A logFWHM=0.13 dex and 
A log ALa=0.4 dex, in order to preserve sufficient statistics 
per box without reducing resolution significantly. The un- 
certainty on the best-fitting parameters induced by a differ- 
ent boxing is negligible with respect to the errors evalued 
through the Monte Carlo procedure. We also verified that 
the results are independent of the adopted division in red- 
shift bins. 

Fig.|3]shows the best- fitting Monte Carlo simulated dis- 
tributions compared to the observed distributions of quasars 
and it is apparent that the fit is good in the allowed region 
of the plane. Table [2] contains the best-fitting values of the 
free parameters, their errors and the corresponding xt va l~ 



■ 

I 



\iz> = 6'4d 

V 

it 



•\ 

f . \ / 

f v 
'.' / \ 

\ - 

'■/': \ 



Kzi>=0.78 

'V:. 



'Viz> = 6'.5d 
\ 

- \ 

■ * v 



\ 



>*z>=0.fl3 
\ i 
\ 

A 

I V 
\ 

\ / 

\ S-vA 

i 



Vz>=6'.6'4' ' ' 1 

\ 
A 
- \ 
" A 

A / 

V 

: : \ 

\ - 
\ 



$z>=088 
\ ! 



A 
: \ 

/ \ 



44 45 46 47 44 45 46 47 44 45 46 47 44 45 46 47 
logXL x (5100 A) [erg/s] 

Figure 4. The 8 panels show the H/3 sample in the FWHM- 
luminosity plane at increasing redshift. Dotted, dashed and dash- 
dotted lines represent the loci of constant monochromatic lumi- 
nosity, constant mass and constant Eddington ratio respectively. 



ues, which indicate that the quality of the fit is comparable 
with that of the Mgn sample. 

The results presented here were obtained under the hy- 
pothesis that the anomalies of the distribution of quasars 
in the Civ FWHM - nuclear luminosity plane are basically 
due to a flattened geometry of the corresponding BLR. In 
order to test this assumption, we re-derived the probabil- 
ity density Pj,fwhm(J, FWHM) convolving Eq. 4 of Paper 
I with the orientation effects predicted under the hypothe- 
sis of a geometrically thin broad line region, in which the 
disc-like component of the gas velocity field accounts for 
the 90 per cent, and assuming that an AGN appears as a 
quasar if it is observed with an inclination angle 8 < 50°. 
Within this assumptions, and adopting the previously de- 
rived values of the free parameters (Tabled the quality of 
the description of the observed distribution of quasars in the 
Civ FWHM - nuclear luminosity plane is good at any line- 
width, indicating that the observed Civ FWHM distribution 
is well consistent with being ascribed to orientation effects. 
However, the best-fitting procedure becomes impracticable 
with the P(,fwhm(^, FWHM) function considered here, be- 
cause it would require to introduce a large number of free 
parameters (e.g. the fraction corresponding to the disc-like 
component of the gas velocity field, the maximum value of 
the inclination angle and their uncertainties). Therefore, the 
reduction of the fitting area to FWHM values greater than 
5000 km/s appears to be the only viable way to proceed 
with present data. 



3.3 Analysis of the H/3 sample 

In Fig. |4]we show the H/3 sample, divided in 8 redshift bins, 
in the FWHM-luminosity plane. We now apply the pro- 
cedure described in Paper I to this sample. Again we dis- 
cretize the distributions in 100 boxes (A logFWHM=0.13 
dex; A log XLx~0A dex). We let i m i n vary as a free param- 
eter in all the redshift bins. 

In Fig. [5] the observed distributions of objects at each 
redshift are compared to the distribution simulated with our 
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Figure 5. The 8 panels show the H/3 sample in the FWHM— luminosity plane at increasing redshift: solid black contour plot (levels: 40, 
110, 200, 300 objects per box, see text) represents the discrete observed distribution of objects. Dotted red contour plot (same levels) shows 
the discrete distribution of a sample of objects simulated with the Monte Carlo method, adopting the assumed P; fwhmCj FWHM) 
probability density (Eq. 4 of Paper I) with the best-fitting parameters (see Table [3). Dotted, dashed and dash-dotted lines represent 
'min, "imax and e max respectively. The rectangle in the first redshift bin represents the boxing adopted to discretize the distributions. 



Table 3. Best-fitting values of minimum luminosity, maximum mass, maximum Eddington ratio and widths of the 
corresponding distributions, with errors and xtt f° r the H/3 sample. In all the redshift bins, the number of degrees of 
freedom is v = 94 (100 data-points and 6 free parameters). 



Bin 


< z > 


#(Objects) 


^min 




^max 


<T m 


6 max 


&e 


X..H 


1st 


0.31 


2689 


43.79±0.01 


0.10±0.01 


9.10±0.15 


0.35±0.05 


-0.35±0.08 


0.23±0.02 


13.0 


2nd 


0.46 


2464 


43.99±0.02 


0.20±0.03 


9.20±0.15 


0.31±0.06 


-0.31±0.07 


0.23±0.03 


13.9 


3rd 


0.56 


2008 


44.10±0.03 


0.24±0.03 


9.25±0.13 


0.35±0.05 


-0.31±0.06 


0.23±0.02 


6.7 


4th 


0.64 


1428 


44.34±0.03 


0.21±0.02 


9.30±0.14 


0.30±0.04 


-0.34±0.07 


0.22±0.02 


9.0 


5th 


0.71 


1256 


44.50±0.03 


0.20±0.01 


9.30±0.10 


0.32±0.03 


-0.35±0.04 


0.22±0.01 


29.6 


6th 


0.78 


834 


44.49±0.02 


0.25±0.01 


9.30±0.08 


0.30±0.02 


-0.35±0.03 


0.22±0.01 


9.5 


7th 


0.83 


392 


44.66±0.02 


0.20±0.02 


9.30±0.09 


0.32±0.03 


-0.21±0.04 


0.22±0.01 


13.7 


8th 


0.88 


266 


44.71±0.03 


0.21±0.02 


9.35±0.11 


0.30±0.03 


-0.17±0.05 


0.22±0.01 


24.2 



Note. Data that come from a best-fitting procedure are displayed in boldface. 



procedure, and Table [2] contains the best-fitting values of the 
free parameters and the corresponding xt values, that indi- 
cate that the quality of the fit is lower than that of the Mgn 
sample. In the following, when comparing the H/3 maximum 
mass estimates with those based on Mgn and Civ data, this 
will be taken into account as the errors derived for the H/3 
best-fitting determinations are somewhat larger. 

Even if there are indications that the BLR clouds where 
the H/3 line is produced are substantially isotropic, as sug- 
gested by consideration on the emission line shape (on this 
point see also the thick disc model in Collin et al. 2006; De- 



car li et al. 2008), we test whether the reduction of the fitted 
region to high FWHM values (see Section [3.2p would allow 
a better agreement between the simulated and the observed 
distribution: It is not the case, confirming that the explana- 
tion of the moderate quality fit of the H/3 data is probably 
uncorrelated with the inclination angles. 

We also test whether the lower quality of the results is 
due to the fact that the H/3 sample is limited to low- z values 
and covers a small redshift window. In order to rule out also 
this possibility, we consider a subsample of the Mgn dataset 
assuming z < 0.91, so that a similar redshift range is covered 
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by both the samples (see Table[T} • The Mgn matched sample 
is split in 8 redshift bins, accordingly to the division defined 
for the H/3 dataset. The application of the maximizing mass 
procedure to this subsample gives \v values smaller than 5 
in each bin, consistently with those obtained for the Mgn 
sample (Table 1 of Paper I) and significantly smaller than 
those derived for the H/3 sample, in the same redshift win- 
dow. Since the objects included in the Mgn matched sample 
coincide with those appearing in the H/3 one, apart for the 
first redshift bin which results less populated, this test in- 
dicates that the anomalies of the fit procedure reside right 
in the H/3-based BH mass estimates (difficulties in spectral 
measurements or exoticism of the BLR geometry). 

Comparing the H/3 to the Mgn FWHM distributions 
for the objects for which both lines are available, we see 
that both the high and the low velocity tails of the H/3 
FWHM distribution are substantially enhanced with respect 
to the Mgn one, although the widths of the two emission 
lines have comparable mean values and correlate reason- 
ably (e.g. McLure & Jarvis 2002, Shen et al. 2008). The 
standard deviation of the H/3 FWHMs (2525 km/s) is con- 
siderably higher than that of the Mgn values (1698 km/s). 
On the other hand, for the same objects, the distribution 
of AL A (3000A) is well consistent with that of AL A (5100A), 
assuming a spectral index a ~ —0.3 (Labita et al. 2008) 
to compare the luminosities at different wavelengths (see 
Fig. [6} ■ This suggests that the anomalies of our results from 
the H/3 sample reside in particular in the odd distribution 
of the H/3 line- widths. 

We then check whether the FWHM values provided by 
Shen et al. (2008), which come from an automatized rou- 
tine, are subject to systematic errors and verified the H/3 
line width measurements for a subsample of randomly se- 
lected objects. In general we find a good agreement, apart 
from some cases in which the FWHM measured by Shen et 
al. (2008) is very high (~ 12000 km/s), while we obtained 
considerably smaller values. This is probably due to diffi- 
culties in the deconvolution of the broad and narrow line 
components, which is a problem peculiar to the H/3 line. 
However, this problem affects only a small number of ob- 
jects, and the dispersion of the H/3 line widths is likely really 
higher than that of the Mgn FWHMs. This would imply that 
the H/3 and Mgn emitting regions are intrinsically different. 
Note that this scenario is hardly consistent with the widely 
accepted assumption that the H/3 and Mgn lines are pro- 
duced in the same region, given that they have similar ion- 
izing potentials (e.g. Verner et al. 1996), unless we suppose 
that the H/3 line is produced in clouds with lower column 
densities with respect to the Mgn line. This would imply 
that the effects of radiation pressure are stronger and that 
the H/3 emitting clouds may be gravitationally unbounded 
(Marconi et al. 2008), accounting for non virial motion ob- 
served from the H/3 emission. 

To explain the larger H/3 FWHM standard deviation 
with respect to the Mgn one, a reasonable possibility is then 
to consider that velocities due to winds, turbulent motion 
or inflowing gas, comparable to the Keplerian ones, may be 
present in the H/3 emitting region (see e.g. Ruff 2008; Mar- 
coni et al. 2008; Gaskell, Goosmann & Klimek 2008 for mod- 
els and observational evidences) . This would also account for 
the mediocre description of the H/3 sample distribution in 
the FWHM-luminosity plane in terms of Eq. 4 of Paper I, 
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Figure 6. The objects common to the Mgn and H/3 sample with 
0.5 < z < 0.7 are plotted in the FWHM-luminosity plane. The 
FWHM and luminosity distributions are also shown (arbitrary 
units) for both Mgn and H/3 data. Vertical thick lines correspond 
to the ±2<r spread in FWHM and horizontal thick lines indi- 
cate the ±2<t dispersion in luminosity. In right panel, thin lines 
trace the FWHM and luminosity dispersion of the Mgn panel, 
assuming a spectral index a ~ —0.3 (Labita et al. 2008) to derive 
logAL A (5100A) from log AL A (3000A). 



as the assumption of the applicability of the virial theorem 
is a fundamental point of our procedure. 



4 QUASAR BH MASS DEPENDENCE ON 
REDSHIFT 

4.1 Combining the H/3, Mgn and Civ samples 

The maximum masses in various redshift bins calculated 
from the three samples are reported in Fig. [7] (upper panel). 
The M B a B ^ i (z) trend appears remarkably smooth, apart from 
modest offsets between the BH mass estimates derived from 
different lines. 

The origin of the two offsets, which are not present in 
the average values per redshift bin of the BH masses from 
Shen et al. (2008), is probably different. In the case of the 
H/3- and Mgll-based mass estimates (hereafter, Mbh ,h/3 and 
AfBH,M g ii), it is principally due to the fact that our points 
are representative of the maximum mass and that the H/3 
and Mgn FWHM distributions are different, in the sense 
that the H/3 FWHM standard deviation is larger than the 
Mgn one (see Section \3. 31 and Fig. [6}. This implies that, for 
objects with large line- widths (and then likely high masses), 
Mbh, h/3 is on average larger than MaH.Mgii- From the data 
by Shen et al. (2008), in fact, we obtain that the difference 
of the BH mass mean values of the objects with both H/3 
and Mgn estimates is negligible (AA/bh = 0.03 dex), while 
it is AA/bh = 0.13 dex in the subsample with both Mbh, h/3 
and Mbh ,Mgii greater than 10 9 Mq. 

On the other hand, the offset between the Mgn- and 
Civ-based BH masses is connected to the choice of the 
FWHM threshold adopted in the analysis of the Civ sample 
(Section l3.2[) . Assuming a higher FWHM limit, the resulting 
Timax estimates would be larger (for instance, it is possible 
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Figure 7. Upper panel: small dots are the the virial BH masses from Shen et al. (2008); the dash-dotted line reports the corresponding 
mean values. Red squares are our estimates of log f rom the H/3 sample (Table [3] column m max ); green triangles are the same 

estimates from the Mgn sample (table 1 of Paper I, column m max ) and blue circles are from the Civ sample (Tableful column m max ). 
Vertical dotted lines indicate the redshift regions in which we derived reliable H/3, Mgn and Civ m max estimates. Lower panel: same 
legend as the upper panel, but the maximum mass estimates from H/3 and Civ were matched to the Mgn estimates (see text). The 
corresponding red, green and blue solid lines are the best linear fit to the data from the H/3, Mgn and Civ samples (Table [4} . The black 
dashed line is the best fit to all the data-points (Eq. 



to cancel the offset with a threshold at 7000 km/s), but the 
trend in redshift remains unchanged, although statistics are 
significantly impoverished. To further test the effects of an 
FWHM cut on the derived redshift dependence of the maxi- 
mum black hole masses, we reproduced the entire procedure 
to the Mgn sample studied in Paper I, reducing the analy- 
sis to objects with Mgn line- widths larger than 5000 km/s. 
We verified that the derived slope of the log M m | x (2) depen- 
dence is consistent with that of Eq. 11 of Paper I within la, 
thus confirming that it is independent of the FWHM cut. 

For these reasons, we apply a rigid shift to the H/3 and 
Civ m max values, imposing that they match the Mgll-based 
estimates in the common redshift regions. The shift amounts 
to -0.06 dex for the H/3 data-points and to +0.16 dex for the 
Civ maximum mass estimates. 



Table 4. Best linear fit to the maximum mass as a function of 
redshift for the H/3 sample, the Mgn sample and the Civ sample. 
The corresponding xl values and the probability of getting worse 
results is also given. 



Sample 


loe" ^ BH 
a 


= az + /3 


v xl 


P(xi 2 > xl) 


H/3 


0.35 ±0.05 


8.97 ±0.04 


6 0.05 


> 99% 


Mgn 


0.34 ±0.02 


8.99 ± 0.03 


5 0.29 


~92% 


Civ 


0.13 ±0.02 


9.45 ± 0.05 


6 0.10 


~99% 
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4.2 A/bh(z) relationship from z ~ 0.1 to z ~ 4 

The shifted values of m max (see previous section) are plotted 
in Fig. [7] (lower panel) as a function of redshift. In Table [4] 
we report the best linear weighted fits from the H/3 sample 
(Section l3.3p . the Mgn sample (Paper I) and the Civ sample 
(Section 13. 2|) . It is worth to notice that the xt of the three 
best-fitting lines are all smaller than 1 and this may indicate 
that the Monte Carlo procedure overestimates the errors on 
the m max values. Rather, the uncertainties that we propose 
are representative of the quality of the description of the 
observed distribution of quasars in the FWHM-luminosity 
plane. 

As the slopes of the log(MBH) - z dependence decrease 
at increasing redshift, i.e. from the H/3 to the Mgll and to 
the Civ sample, we choose to describe the overall trend of 
the maximum quasar BH mass estimates with a power law. 
The best fit is: 



M. 



oc (1 + z) 1 



(2) 



We stress that the importance of Eq. [2] is limited to the 
resulting trend in redshift rather than to the normalization 
factor, which comes from a number of assumptions on the 
definition of 'maximum mass'. 

The overall best-fitting function is practically indistin- 
guishable from Eq. 11 of Paper I in the Mgll redshift range 
(see Figure lower panel). We thus confirm and extend the 
results of Paper I on a wider redshift range. 

Assuming that the shapes of Mbh and Eddington ratio 
distributions do not change appreciably with redshift, which 
is suggested by the fact that the value of a m are practically 
independent of z for each virial estimator (see Tables [3] [2] 
and table 1 of Paper I), Eq.[2]also describes the trend in red- 
shift of the mean quasar BH mass. We conclude that quasar 
samples at higher z are increasingly dominated by higher 
mass BHs, i.e. the lower is the mass the longer is the cosmic 
age at which quasars exhibit their activity. More precisely, 
we find evidence that the most massive black holes stop nu- 
clear activity before less massive ones, since we describe the 
upper envelope of the BH mass distribution. In any case, this 
is indicative of anti-hierarchical 'cosmic downsizing', which 
refers to the changing epoch of the peak of quasar activity 
for high and low mass black holes. 



4.3 Comparison with previous results 

In Fig. [8] we present a global comparison between our re- 
sults and the state-of-art knowledge on the redshift evolu- 
tion of the quasar black hole mass. Our findings are consis- 
tent within the la uncertainty with those by other authors 
(e.g. Fine et al. 2006 and Decarli et al. in preparation, see 
Section 4.3 of Paper I for a detailed discussion). Our results 
are in agreement with recent theoretical predictions by e.g. 
Merloni (2004) and Merloni, Rudnick & Di Matteo (2006), 
who studied the redshift evolution of the BH mass func- 
tion and computed the average mass of supermassive active 
black holes as a function of z through models accounting 
for radiatively efficient and inefficient BH accretion. There 
is also consistency with models by e.g. Shankar, Weinberg 
& Miralda-Escude (2009), who studied AGN duty-cycles 
through BH accretion models to predict the z-dependence 



of the minimum BH mass at which supermassive black holes 
start to exhibit nuclear activity. The results by Marconi 
et al. (2004), who modeled with similar methods the z- 
dependence of the BH mass of active BHs that accreted a 
fixed fraction of their total mass (i.e., which entered their 
AGN phase) , are also qualitatively consistent with our find- 
ings. 

Fig. [8] also shows the Mbh(z) trend corresponding to a 
minimum detected flux, which is representative of the effects 
of a Malmquist-type bias on the observed average BH mass 
redshift dependence. In detail, we considered a survey flux 
limit that roughly corresponds to that of the SDSS quasar 
selection algorithm (i < 19.1 at z < 3 and i < 20.2 at z > 3, 
see Section 13.21 and references therein) . From the luminos- 
ity dependence of the virial BH mass estimator (see Sec- 
tion 2) we estimate the corresponding mass limit at a given 
FWHM as a function of redshift. It is apparent from Fig. [8] 
that the derived BH mass trend traces the observed aver- 
age Mbh z-dependence of the SDSS quasars at any redshift, 
thus confirming that it is dominated by the Malmquist bias. 
On the other hand, our results deviate substantially from 
the predicted effects of the bias, particularly at low redshift. 



5 SUMMARY AND DISCUSSION 

We presented the application of a procedure which allows 
the determination of the z-dependence of the BH mass to a 
large sample of quasars for which line-width and continuum 
measures are available. 

Combining the results from the H/3, Mgll and Civ sam- 
ples of SDSS quasars, we found that the maximum mass of 
the active BH populations evolves as (1 + 2 ) 1 ' 64±0 04 from 
z ~ 0.1 to z ~ 4. It is noticeable that the dependence on 
redshift of the typical mass of active black hole populations 
is substantially unchanged from ^ of the actual age of the 
Universe to now, covering the region of maximum quasar 
activity (z ~ 2 — 3). This is a clear manifestation of cosmic 
downsizing: the most massive black holes were more actively 
accreting early in the Universe and with time progressively 
less massive ones dominate the population of AGN. 

A comparison between the three emission lines stud- 
ied in the work suggests some considerations on the corre- 
sponding BLR properties. Consistently with other indepen- 
dent indicators of the BLR geometry (considerations on the 
broad lines shape and from the geometrical factor /), we 
suggest that the Mgll BLR is essentially isotropic (or rather 
geometrically thick) and the Doppler velocities are practi- 
cally purely rotational, as the assumptions produce a high 
quality description of the data distribution in the FWHM- 
luminosity plane. Consequently the Mgll line width is a good 
indicator of the virial velocity from which it is possible to 
infer the BH mass. On the other hand, there is increasing 
evidence that the Civ line emitting region is disc-like and 
hence this line is representative of the virial velocity only 
if the BLR is observed edge-on, i.e. when the Doppler ef- 
fect (which acts only on the velocity component parallel to 
the line of sight) is maximum. The BH mass determinations 
from Civ are then subject to an intrinsic error due to the 
uncertainty on the inclination angle and then on the virial 
velocity of the clouds. The H/3 line emitting region of the 
BLR is found to be closer to the isotropic case than Civ 
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Figure 8. Thin solid line is the M^J x (z) dependence studied in 
this work (Eq. [2}- Dot-long dashed line is the best fit of Mbh( z ) 
by Fine et al. (2006) for a sample of 2dF quasars with luminosity 
around L*(z) at each redshift. Long dashed line is the best fit 
of Mbh as a function of redshift for a sample of ~ 100 quasars 
with high quality spectroscopic observations (Decarli et al. 2009, 
in preparation). Shaded areas represent the uncertainties of the 
relationships. Dotted line shows the average Mbh as a function 
of redshift of the SDSS DR5 quasars (Shen et al. 2008) and the 
corresponding shaded area illustrates the error on the mean values 
per redshift bin. Green thick line represents the Mbh (2) trend 
corresponding to the minimum flux of the SDSS quasar selection 
algorithm (see text). Red lines are theoretical predictions from BH 
accretion models: dot-short dashed line is the average Mbh (2) by 
Merloni et al. (2004), short-long dashed line represents Mbh(z) 
of BHs that accreted 5 per cent of their total mass (Marconi et 
al. 2004); short dashed line is the minimum BH mass required to 
exhibit nuclear activity as a function of z (Shankar et al. 2009). 



one, but we suggest that for some objects the H/3 clouds are 
not at the equilibrium and are dominated by turbulent mo- 
tion. This indicates that the H/3 emission line is a mediocre 
estimator of the BLR virial velocity and that the BH mass 
estimates from H/3 are subject to an uncertainty which is 
difficult to constrain. 

In conclusion, we suggest that an inflated disc broad- 
line region, in which the Carbon line is emitted in a flat inner 
disc while Mgn is produced in a geometrically thick outer re- 
gion, can account for the observed differences in the FWHM 
distributions of the two emission lines. The H/3 broad emis- 
sion line would be produced in clouds at the same distance 
from the central source as the Mgn ones, but the role of ra- 
diation pressure may be more relevant, accounting for non 
virial motion of the gas. We then propose that the Mgn line 
and the corresponding 3000 A monochromatic luminosity 
are preferable indicators of the BLR virial velocity and size 
respectively, and hence of the BH mass through Eq. [T] 
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